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A linked cluster expansion for the calculation of ground state observables of 
' complex nuclei with realistic interactions has been used to calculate the ground 

state energy, density and momentum distribution of ^^O and ''"Ca. Using the 
same cluster expansion and the wave function and correlation parameters ob- 
tained from the energy calculation, we have evaluated the semi inclusive reaction 
A{e, e'p)X taking final state interaction (FSI) into account by a Glauber type 
\ approach; the comparison between the distorted and undistorted momentum 

distributions provides an estimate of the transparency of the nuclear medium 
to the propagation of the hit proton. The effect of color transparency is also 
■ included by considering the Finite Formation Time (FFT) that the hit hadron 

' needs to reach its asymptotic physical state. 

m ' 
o 

1 Introduction 

The exclusive, A{e, e'p)B, and semi-inclusive, A{e, e'p)X, electro-disintegration 
y . of nuclei represent a powerful tool to investigate various aspects of nuclear struc- 

^ I ture (e.g. single-particle motion and mean field effects, nuclear correlations), as 

well as QCD motivated effects (e.g. color transparency). The accuracy of recent 
and forthcoming experimental data require realistic theoretical calculations to 
^ . be performed, based as much as possible on a description of the nucleus stem- 

$H ' ming from first principle calculations, which means that the nuclear wave func- 

tions which appear in the calculation of various matrix elements, used either to 
predict or interpret the experimental data, should in principle result from many 
body calculations and realistic interactions. The problem has been solved in the 
case of few body system, for which realistic wave functions are currently being 
used in the interpretation of electro-disintegration processes, but this problem 
still needs a solution in case of complex nuclei. As a matter of fact, calculations 
of ground state observables for complex nuclei represent still a hard task, and 
even in those cases when approximate many body calculations for the ground 
state energy can be performed, the structure of the wave function is so compli- 
cated that its use for calculations of matrix elements of operators, different from 
the potential and kinetic energies ones, is very involved. For such a reason, a 

^Presented by M. Alvioli at the 2nd International Conference on Nuclear and Par- 
ticle Physics with CEBAF at JLab, to appear in the Proceedings. 
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simpler, but still realistic, method which would allow one to calculate various 
kinds of matrix elements with nuclear wave functions which correctly incorpo- 
rate the most relevant features of a realistic wave function, in particular its 
correlation structure resulting from the main features of modern two-nucleon 
interactions, would be extremely useful. Cluster expansion techniques, when 
the expectation value of different operators can be calculated to a certain order, 
may represent a valid and practicable alternative to the full complex "exact" 
solution of the many body problem. It is the aim of the present paper to il- 
lustrate a cluster expansion approach to the calculation of ground state energy 
properties (energy, density and momentum distribution) and various types of 
electro-disintegration processes. 



2 Cluster expansion and the nuclear wave func- 
tion 

In our linked-cluster expansion approach, the expectation value of a certain 
operator O 

is evaluated with correlated wave functions of the following " classical" form 

- F(ri,...,r^)$A(ri,...,r^), (2) 

where is a mean field (Slater determinant) wave function, and F a sym- 
metrized (by the symmetrization operator S) correlation operator which gener- 
ates correlations into the mean field wave function; it has the following general 
form 

A 

with 

= E /^'^(^'^■) z^'H^^^) = /^'^(^^^■) (4) 

p 

where the operators O'-P-' are the same which appear in the two-nuclcon inter- 
action, having the form ( e.g. in case of a VS-type interaction) 

Off''' = [1, ^^ ■ o,, S,j, (L • S),,] ® [1, n ■ T,] (5) 

The central parts /^P-'(ry )'s of the correlation function /'■P-', reflect the radial 
behaviour of the various components and their actual form is determined either 
by the minimization of the ground state energy, or by other criteria. 

The cluster expansion of Eq^ is carried out in terms of the quantity fjij = 
ffj — 1, whose integral plays the role of a small expansion parameter; by ex- 
panding the numerator and the denominator the terms 0„ of the same order n 
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in r]ij, are collected obtaining (O) = On + Oi + O2 + with 

ij ij 

02 ^ {J2 fi^^ ^) - (E ^) ^'^) + 

ij<kl ij ij 

+ C'o ( E '^'^^ ) 

\ ij<kl 

where ([...]) = ($a |[---]| ^a)- From now on, our approach will consist in obtain- 
ing the parameters characterizing the correlation functions and the mean-field 
single-particle wave function which correspond to an acceptable value of the 
ground state energy, we will then use the obtained wave function ^ a to cal- 
culate the transition matrix elements entering in the theoretical description of 
electro-disintegration processes using the same cluster expansion employed to 
calculate the energy. We have calculated the ground state energy of ^^O and 
^°Ca using the Argonne V%' P potential and adopting, as in Ref. [2], the so 
called /e approximation consisting in considering only the first six components 
of Eq. [S] The expectation value of the many body non relativistic Hamiltonian 
of the nucleus was obtained by calculating the average values of the kinetic and 
potential energies, i.e. 

(T) = -|^/ dkfc2n(A:), (7) 
where n(fc) is the nucleon momentum distribution (fc = |k|), 

n{k) = / rfridrle-^'' (-^-"^i)p(i)(ri,ri), (8) 

and 

kV) - ^E(^^^-)= ^^^^E / dridr2«(^)(ri2)pg;(ri,r2). (9) 

The calculations have been performed by cluster expanding the expectation 
value of the non diagonal one-body, /5^^\ and diagonal two-body, p^^^(v\,V2) 
density matrix operators. The six correlation functions /^^H^ij) have been 
borrowed from Ref. [2], whereas Harmonic Oscillator (HO) and Saxon- Woods 
(SW) spwf's have been used to describe the mean field. As in Ref. [5], we 
found that the charge densities corresponding to the minimum of the energy, 
appreciably disagree with the corresponding experimental quantities, therefore, 
in view of the mild dependence of the energy around the minimum upon the 
mean-field parameters, following Ref. |2] we have changed the latter to obtain 
agreement between theoretical and experimental charge densities. The results 
for the charge densities and momentum distributions, which are shown in Figs. 
n]and[21 deserve the following comments: 
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1. the agreement between our cluster expansion and FHC/SOC result of Ref. 

is very good; 

2. both approaches predict momentum distributions which do not apprecia- 
bly differ from the ones obtained in Ref. 3 , where the Variational Monte 
Carlo method and the AV^18 interaction have been used; 

3. the high momentum part of n{k) is almost entirely exhausted by non- 
central, long-range correlations, with the central, short-range, Jastrow cor- 
relations under-predicting the high momentum part of n{k) by about one 
order of magnitude; 

4. the dominant non-central correlations are the isospin, /4 — f^^\rij)Ti ■ tj, 
and isospin-tensor, /g — f'^^\rij)Ti ■ TjSij, correlations; 




Figure 1: The charge densities of ^^O calculated using the cluster expansion 
lO with Harmonic Oscillator (HO) (left panel) and Saxon- Woods (SW) (right 
panel) single-particle wave functions {spwfj. Dashed line: mean- field wave func- 
tions used in the cluster expansion calculations; full line: results of the cluster 
expansion; thick full line: experimental data j^. The HO and the SW spwf's 
parameters corresponding to the full line have been chosen such as to pro- 
vide a good description of the density without appreciably changing the min- 
imum value of the ground state energy. The normalization of the density is 
477 / p(r) dr = Z, Z being the number of protons. 

3 The final state interaction in A{e, e'p)X reac- 
tions off complex nuclei: Glauber approach 

Using the results obtained in the previous Section, we have calculated the 
semi-inclusive A(e,e'p)X process in which an electron with 4-momentum ki = 
{ki,zei}, is scattered off a nucleus with 4-momentum Pa = {0, zM^i} to a 
state k2 = {k2,ie2} and is detected in coincidence with a proton p with 4- 
momentum p = {p,iEp}; the final {A — 1) nuclear system with 4-momentum 



4 




Figure 2: The momentum distributions of ^^O calculated with the cluster ex- 
pansion (|SJ) and Harmonic Oscillator (HO) (upper panel) and Saxon- Woods 
(SW) (lower panel) spwf's as in Fig.l. Dashed line: mean field approximation; 
thick full line: fi approximation (only central correlations: (/ = /^^■'); dotted 
line: /a approximation (/ — f'^^^ + + /'^•'); full line: fg approximation 
if = ELi /^^^)- The asterisks and the fuU squares are the results from Ref. 
121 and 0, respectively. The values of the kinetic energy calculated by Eq. [7| 
are as follows: < T > = 230.28 MeV (dashes), 287.23 MeV (thick full), 435.94 
MeV (dots), 458.29 MeV (full) for H.O. spwf's, and < T > = 244.94 MeV 
(dashes) 306.98 MeV (thick full), 469.83 (full), 494.47 MeV (dots) for S.W. 
spwf's. The normalization of n(k) is 4tt J n(k) fc^ dk — 1. 



Px = {PxjiEx} is undetected. The cross section for the exclusive process 
A(e,e'p)B can be written as follows 

'^'^ = KaepPD(Em,Pm) (10) 



dQ^dvdp 
d factor, o 

and = |qp — v'^ the four momentum transfer. The quantity PD(Em,Pm) 



where if is a kinematical factor, a^p the off-shell electron-nucleon cross section, 
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is the distorted nucleon spectral function which depends upon the observable 
missing momentum p,„ = q — p (Pm = k when the FSI is absent) and missing 
energy E,n — v — Tp~ Ta-i- In the semi- inclusive A{e, e'p)X process, the cross 
section (|10|) is integrated over the missing energy Em , at fixed value of Pm and 
becomes directly proportional to the distorted momentum distribution 

noiPm) = (2^)-3y"e^P'"("-i-"-i)pc(ri,r;)dridr; (11) 

where 

^-(^I'^i)- 1^^) (12) 

is the distorted one-body mixed density matrix, S the S-matrix describing FSI, 
and the primed quantities have to be evaluated at with i — 1,...,^. The 
integral of n/)(p,„) gives the nuclear transparency T 



T . 



J nD{Pm)dp 
/ n{k)dk 



= J pDir)dr = 1 + AT (13) 



where Pd{j^) — Pd{^i — = r) and AT originates from the FSI. In Ref. |1] Eq. 
II II has been evaluated using a Glauber representation for the scattering matrix 
5*, viz 

A A 

S ^ Sain . . . r^) = n G(ri, r,) 0{z, - z^)r{h, - b,)] (14) 

where and Zj are the transverse and the longitudinal components of the 
nucleon coordinate = {hj,Zj), r{h) the Glauber profile function for elastic 
proton nucleon scattering, and the function 9{zj — zi) takes care of the fact 
that the struck proton "1" propagates along a straight-path trajectory so that it 
interacts with nucleon "j" only if z^- > zi. The same cluster expansion described 
in Section II has been used taking Glauber rescattering exactly into account at 
the given order n, and using the approximation I^'a-sP = lis Using the 

mean-field and correlation parameters obtained from the energy calculation, 
we have obtained the distorted nucleon momentum distributions nD{pm) = 
nD{pm, 0), where 9 is the angle between q and Pm', the results for ^^O and "^^Ca 
are presented in Fig. |2| 



4 Finite formation time effects 

Recently [B] the effects of color transparency in quasi-elastic lepton scattering off 
nuclei, have been introduced by explicitly considering the finite formation time 
(FFT) that the hit hadron needs to evolve to its asymptotic physical state. It has 
been shown that at the values of the Bjorken scaling variable x — Q'^ /2mv ~ 1, 
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Figure 3: The distorted momentum distribution, no{Pm) = nD[pm,6) [6 = 
c^j„J, of ^^O and ^'^Ca, obtained by Eq. ^2 with correlated wave functions, 
Harmonic Oscillator spwf's and the Glauber S matrix p4|l . The value of the 
integrated nuclear transparency H13|l for ^^O is 0.5. 
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Figure 4: Left panel: the Q^- and Pm-dependence of FFT effects on the distorted 
momentum distribution of ^^O calculated by Eqs. 1111 and [HI at 6 — 90°. Right 
panel: the ratio of the distorted to the undistorted momentum distributions, i.e. 
the non- integrated nuclear transparency for various values of Pm and 9 — 90°. 
The horizontal lines represent the Q^-independent Glauber results. 



FFT effects can be treated in a simple way, i.e. by replacing the Glauber 
operator (Eq. I14|l with 

A 

SFFT{ru...,rA)^ n(l- J(zi-z,)r(bi-b,)), (15) 

where 

J{z) = 6l(z)e"^ (16) 

m being the nucleon mass and = m* ^ — m? is a parameter describing the 
average excitation energy of the ejectile. It can be seen that at sufficiently high 
values of Q^, J — > 1 and the FSI vanishes. The effects of FFT on the distorted 
momentum distribution is shown in Fig. 01 
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5 Summary and conclusions 



We have obtained fully correlated wave functions by calculating the average 
value of the nuclear Hamiltonian by means of a linked cluster expansion and 
using realistic two-nuclcon interactions. The wave functions have been used to 
obtain the ground state density and momentum distribution. By introducing 
FSI effect by a Glauber-type approach, the distorted momentum distributions 
appearing in the semi-inclusive A(e, e'p)X processes have been calculated. By 
such a procedure, a consistent treatment of initial-state correlations and final- 
state interactions has been achieved. Color transparency effects have also been 
investigated by the Finite Formation Time approach. Comparison with avail- 
able experimental data are in progress and will be reported elsewhere. 

This work is part of the research activity of M. A. performed, under the 
supervision of C. d. A., for the fulfillment of the PhD title. H. M. thanks the 
Department of Physics, University of Perugia, and INFN, Sezione di Perugia, for 
hospitality and support. Partial support by the Italian Ministero dellTstruzione, 
Universita e Ricerca (MIUR), through the funds COFINOl, is acknowledged. We 
are indebted to A. Fabrocini for providing the correlation functions obtained in 
Ref. 12 . 
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